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Abstract. Hydrostatic equilibrium and energy conservation determine the 
conditions in the gravitationally stabilized solar fusion reactor. We assume a 
matter density distribution varying non-linearly through the central region of 
the Sun. The analytic solutions of the differential equations of mass conserva- 
tion, hydrostatic equilibrium, and energy conservation, together with the equa- 
tion of state of the perfect gas and a nuclear energy generation rate e = eop n T m , 
are given in terms of Gauss' hypergeometric function. This model for the struc- 
ture of the Sun gives the run of density, mass, pressure, temperature, and nuclear 
energy generation through the central region of the Sun. Because of the assump- 
tion of a matter density distribution, the conditions of hydrostatic equilibrium 
and energy conservation are separated from the mode of energy transport in the 
Sun. 



1 Hydrostatic Equilibrium 

In the following we are concerned with the hydrostatic equilibrium of the purely 
gaseous spherical central region of the Sun generating energy by nuclear re- 
actions at a certain rate (Chandrasekhar, 1939/1957; Stein, 1966). For this 
gaseous sphere we assume that the matter density varies non-linearly from the 
center outward, depending on two parameters 5 and 7, 



where x denotes the dimcnsionless distance variable, x = r/R Q , < x < 1, R Q 
is the solar radius, S > 0, 7 > and 7 is kept a positive integer in the following 
considerations. The choice of the density distribution in (1) and (2) reveals 
immediately that p(x = 0) = p c is the central density of the configuration 
and p(x = 1) = is a boundary condition for hydrostatic equilibrium of the 
gaseous configuration. For the range < x < 0.3 the density distribution in 
(1) and (2) can be fit numerically to computed data for solar models by chosing 
5 = 1.28 and 7 = 10 (Haubold and Mathai 1994). For these values of S and 
7 the function fo(x) in (2) is shown in Figure 1. The choice of restricting 
x to x < 0.3 is justified by looking at a Standard Solar Model (Bahcall and 
Pinsonneault, 1992) which shows that x < 0.3 comprises what is considered to 
be the gravitationally stabilized solar fusion reactor. More precisely, 95% of the 
solar luminosity is produced within the region x < 0.2(M < 0.3M Q ). The half- 
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peak value for the matter density occurs at x = 0.1 and the half-peak value for 
the temperature occurs at x = 0.25. The region x < 0.3 is also the place where 
the solar neutrino fluxes are generated. As we are concerned with a spherically 
symmetrical distribution of matter, the mass M(x) within the radius x having 
the density distribution given in (1) and (2) is 

M(x) = M f M (x), (3) 
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7! do 
where M© denotes the solar mass and 2-F1 (•) is Gauss' hypergeometric function 
(see, for example, Mathai, 1993). Equations (3) and (4) are satisfying the 
boundary condition M(x = 0) = and determine the central value p c of the 
matter density through the boundary condition M(x = 1) = M©, where p c 
depends then only on 5 and 7 of the chosen density distribution in (1) and (2). 
The function /m(#) in (4) is shown in Figure 2, using Mathematica (Wolfram, 
1991). 

For hydrostatic equilibrium of the gaseous configuration the internal pressure 
needs to balance the gravitational attraction. The pressure distribution follows 
by integration of the respective differential equation for hydrostatic equilibrium, 
making use of the density distribution in (1) and the mass distribution in (3), 
that is 
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where G is Newton's constant and 2-^1 (•) denotes again Gauss' hypcrgeometric 
function (Mathai, 1993). 

The Pochhammer symbol (|+m + l) 7 = r(| + m + l + 7)/r(| + m + l) often 
appears in series expansions for hypergeometric functions. Equations (5) and 
(6) give the value of the pressure P c at the centre of the gaseous configuration 
and satisfy the condition P(x = 1) = 0. The graph of the function fp(x) in (6) 
is shown in Figure 3, using Mathematica (Wolfram, 1991). 



2 Equation of State 

It should be noted that P(x) in (5) denotes the total pressure of the gaseous 
configuration, that is the sum of the gas kinetic pressure and the radiation pres- 
sure (according to Stefan-Boltzmann's law)(Chandrasckhar, 1939/1957; Stein, 
1966). However, the radiation pressure, although the ratio of radiation pressure 
to gas pressure increases towards the center of the Sun, remains negligibly small 
in comparison to the gas kinetic pressure. Thus, Equation (5) can be consid- 
ered to represent the run of the gas pressure through the configuration under 
consideration. Further, the matter density is so low that at the temperatures 
involved the material follows the equation of state of the perfect gas. Therefore, 



the temperature distribution throughout the gaseous configuration is given by 

(7) 
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where A: is the Boltzmann constant, Na Avogadro's number, /u, the mean molec- 
ular weight, and 2-^1 (•) Gauss' hypergeometric function (Mathai, 1993). Equa- 
tions (7) and (8) reveal the central temperature for T(x = 0) = T c and satisfy 
the boundary condition T(x = 1) = 0. Since the gas in the central region of the 
Sun can be treated as completely ionised, the mean molecular weight \i is given 
by [i — (2X + jY + where X, Y, Z are relative abundances by mass 

of hydrogen, helium, and heavy elements, respectively, and X + Y + Z = 1. 
The run of the function fr(x) in (8) is shown in Figure 4, using Mathematica 
(Wolfram, 1991). 



3 Nuclear Energy Generation Rate 

Hydrostatic equilibrium and energy conservation are determining the physical 
conditions in the central part of the Sun. In the preceeding Sections, the run of 
density, mass, pressure, and temperature have been given for a gaseous config- 
uration in hydrostatic equilibrium based on the equation of state of the perfect 
gas (Equations (l)-(8)). In the following, a representation for the nuclear energy 



generation rate, 

e(p,T) = e p n T m , (9) 

will be sought which takes into account the above given distributions of density 
and temperature and which can be used to integrate the differential equation of 
energy conservation throughout the gaseous configuration considered in Sections 
1 and 2. In Equation (9), n denotes the density exponent, m the temperature 
exponent, and e a positive constant determined by the specific reactions for 
the generation of nuclear energy. Using the equation of state of the perfect gas, 
Equation (9) can be rewritten more conveniently, 
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where we note that < P{x)/P c < 1, and P(x) is given in (5). Subsequently, 
we can write P r as follows: 
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Taking the ratio of the pressure P(x) at the location x to the central value of 



it one can write 
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Note that h(x) is a polynomial of degree 2j in x s . Denoting the polynomial 
h(x) by 

h(x) = a + ai [x s ] + a 2 [x 5 } 2 + . . . + a 27 [x s }^ , (15) 
we obtain for (13) with a view to (10), 
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where we can expand ^(ai)] 9 by using a multinomial expansion. That is 
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where n + ni + . . . + n 27 = q. Note also that since < x 5 < 1 we have for 
a; < 1, in Equation (10), 
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For the nuclear energy generation rate in (10) we obtain finally 
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where 
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with no + ni + . . . + n2- t = q. 

Equation (19) reflects the analytic representation of the nuclear energy genera- 
tion rate in (9) taking into account the run of physical quantities given for the 
gaseous configuration in hydrostatic equilibrium considered in Sections 1 and 2. 
In deriving the explicit dependence of e on the density and temperature we have 
exercised particular care, because the rate of nuclear energy production is very 
highly temperature sensitive. Small changes in the temperature in the central 
part of the Sun are adequate to balance large differences in luminosity. 

4 Total Nuclear Energy Generation 

The total net rate of nuclear energy generation is equal to the luminosity of the 
Sun, that means the generation of energy by nuclear reactions in the central 
part of the Sun has to continually replenish that energy radiated away at the 
surface. If L{x) denotes the outflow of energy across the spherical surface at 
distance x from the center, then in equilibrium the average energy production 
at distance x is 



where p(x) and e(x) are given in (1) and (19), respectively (Chandrasekhar, 
1939/1957; Stein, 1966). Collecting all factors containing the relative distance 
variable x, the integral to be evaluated is the following, denoting it by g(x), 
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by setting a;" 5 = v. Equation (22) represents an incomplete beta function which 
can be written as a series or in terms of a hypergeometric function. Let s* = 
s + |(3 + 2q) + m + 2n 2 + . . . + (27)n 27 , then we have 

i rx« 



5 (x) = PcR%\J q dv(l-v)->v 



SM 



(-l)i k 
l\ s* + l 

1=0 



1 



= PcR 3 j^[x 5 f 2Fi(-j,s*;s* + l:x s ), (23) 
where 2^1 (•) is Gauss' hypergeometric function (Mathai, 1993). Hence 

L{x) = 4 7 reopr 1 T c m i?|i/l[^] s * 2 F l {- 1 , s*; s* + l;x s ), (24) 
where f is defined in (20). The luminosity, L(x = 1) = L®, is given by 

Lq = 47re p" +1 T" l i? 3 i / — — ^- -, (25) 
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using the relation 2 F 1 (a, b; c; 1) = [T(c)T(c - a - b)]/[T(c - a)T(c - 6)] (Mathai, 
1993). 
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